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ABSTRACT 

The third Progress In Electromagnetics Research Symposium (PIERS) was held 
July 12-16, 1993, at Jet Propulsion Laboratory, California Institute of Technology, 
Pasadena, California. More than 800 presentations were made, and those abstracts are 
included in this publication. 
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SESSION 1A1 

Monday, July 12, AM 0820-1140, Beckman Institute Auditorium 
Organizer: D. Lesselier 

INVERSE PROBLEMS OF LOW-FREQUENCY ELECTROMAGNETISM 

Co-Chairpersons: D. Lesselier, CNRS-ESE, Gif-sur-Yvette, France, and 
H. Ermert, Ruhr— University Bochum, Bochum, FRG 

0820 Inverse Transformation From a Diffusive to a Propagative Field: Analytical 
and Numerical Solutions 

D. Gibert and B. Tournerie, Laboratoire de Geophysique Interne, 

Geosciences Rennes, Universite de Rennes, France 

J. Virieux, Institut de Geodynamique, CNRS, Valbonne, France 8 ~ 

0840 Optimization of Power Devices: A Review 

S. R. H. Hoole, Harvey Mudd College, Claremont, CA, USA......... 9 

0900 A Broadband Holographic Imaging Algorithm for Eddy Current NDT 

H. P. Maier and H. Ermert, Ruhr — University Bochum, Bochum, FRG 

L. von Bemus, Siemens AG, Power Generation Group, Erlangen, FRG 10 

0920 Eddy Current Imaging of Defects in a Conductive Half-Space: Diffraction 
Tomographic and Generalized Solutions in a Wavefieid Framework 
R. de Oliveira-Bohbot, D. Lesselier, B. Duchene, M. Nikolova, and 

A. Mohammad-Djafari, CNRS-ESE, Gif-sur-Yvette, France......... 11 

0940 Considerations on the Determination of the Spatially Varying Conductivity a 

nd Permeability of Metals from Eddy-Current Data 

J. H. Rose, Center for NDE, Iowa State University, Ames, IA, USA — .12 

1000 Coffee Break 

1020 Application of Markov Models for Inverse Problems in Eddy Current 

Nondestructive Evaluation 

L. Udpa and W. Lord, Iowa State University, Ames, IA, USA 13 

1040 Parabolic Electromagnetic NDE Phenomena 

W. Lord, S. Ross, L. Udpa, and M. Lusk, Iowa State University, Ames, IA, USA 14 

1100 An Introduction to Low Frequency (ULF/ELF/VLF) Radio Polarimetry and 

Some Applications to the Interpretation of Natural and Cultural Emission Signatures 
J. Y. Dea and P. M. Hansen, NCCOSC, NRaD, San Diego, CA, USA 

W.-M. Boerner, University of Illinois at Chicago, Chicago, IL, USA 15 

1120 Low Frequency Inverse Scattering 

P.-L. Shen and S.-M. Lin, Northwestern Polytechnical University, Xian, 

Shaanxi, PRC ...16 ' 
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Inverse transformation from a diffusive to a propagative field: 
analytical and numerical solutions. 

Dominique GIBERT*, Jean VIRIEUX 1 ', Benoit TOURNERIE*. 

* Laboratoire de Gdophysique Interne, Geosciences Rennes, University de Rennes 1, Campus de Beaulieu, 35042 Rennes 
Cedex, France. 

t Institut de G£odynamique, CNRS, Avenue A. Einstein, 06560 Valbonne, France. 


Imaging the Earth interior is a major task in geophysics which involves many different techniques 
depending on both the nature and the depth of the targets to detect. For instance, long-period surface 
waves are an efficient tool to probe the upper mantle down to 700 km, while artificially created elastic 
waves constitute the preferred mean to image the first few kilometers of the subsurface in the context 
of petroleum exploration. For high-resolution imaging of the subsurface, the main operational techni- 
ques are seismic (elastic waves) and georadar soundings (electromagnetic waves). 

Unfortunately, while the seismic methods can be adapted to any penetration depth, ohmic dissipa- 
tion strongly limits the georadar penetrating power about 100 meters. Greater depths can be reached 
by using low-frequency electromagnetic fields which do not propagate but rather diffuse in the 
ground. This makes low-frequency electromagnetic probing a notoriously low-resolution imaging 
method. However mathematical links exist between the telegrapher equation, the wave equation, and 
the diffusion equation. For example, the transformation from a propagative to a diffusive field invol- 
ves a Fredholm equation of the first kind 1 . The purpose of our study is to document the inverse of this 
transformation. 

An exact solution for the inverse problem is obtained by performing a spectral decomposition of 
the integral operator 2 . This formal solution gives us a first glance at the ill-posedness of the transfor- 
mation; in particular, the influence of noise is precisely evaluated. If formulated in the context of sharp 
reflector imaging, the inverse transformation benefits from strong a-priori information under the form 
of sparsity constraints. In this framework, the inverse problem takes a form similar to the well-known 
Prony model of spectral analysis. However, the standard numerical solutions of the Prony analysis 
strongly fail when faced with the low signal-to-noise ratio of our data. This comes from the fact that 
these solutions are linearized approximates. We present alternative numerical solutions designed to 
explicitely account for the strong non-linearities of the inverse problem. Although successful! in the 
reconstruction of the reflectors, the simulated algorithm 3 suffers slow convergence rates. A second 
approach uses a coupling of both the simulated annealing and the simplex algorithms, and allows for 
a better account of the continuum nature of the solution space. Finally, we present our first results con- 
cerning the neural network algorithm. All these solutions are tested in the context of two-dimensional 
magnetotelluric soundings with noisy synthetic data. 


1. Lee K.H., G. Liu, and H.F. Morrison, Geophysics, 154, 1180-1192, 1989. 

2. McWhirter J.G., and E.R. Pike, J. Phys., All, 1729-1745, 1978. 

3. Gibert D., and J. Virieux, J. Geophys. Res., B96, 8057-8067, 1991. 
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Optimization of Power Devices: A Review 

S. Ratnajeevan H. Hoole 

Department of Engineering, Harvey Mudd College, Claremont, CA 91711, USA 


Optimization in designing electromagnetic power devices is now increasingly better 
understood. As opposed to classical circuit models of magnetic circuits, today, gradient, stochastic 
and search techniques for mathematical optimization have been proposed and are used. These 
techniques, while being expensive, are exact. More recently, artificial neural networks have been 
suggested, but they work best only if the data set of parameter-set, performance pairs for training 
the network is close to the optimal solution we seek. 

In this paper, the state of the art is reviewed and it is shown how all these methods may be 
used in concert to increase efficiency. The circuit model is used to generate an approximate inverse 
solution. Then direct finite element solutions are used to generate the required training set and this 
is used with the neural network to get a better solution. This solution is finally used as a starting 
point for the gradient optimization scheme which converges quickly because the starting point is 
close to the actual solution. 
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A Broadband Holographic Imaging Algorithm for Eddy Current NDT 
Hans Peter Maier, Helmut Ermert, and Ludwig von Bernus t 


Ruhr-University Bochum, 
Department of Electrical Engineering, 
Univ.-Str. 150 / Geb. IC 
4630 Bochum, FRG, 


t Siemens AG, Power Generation Group 
KWU-S411 

Instrumentation & Testing NDE- Laboratory 
8520 Erlangen, FRG 


The use of eddy currents is a widespread method for nondestructive testing of electrically conductive 
materials. Applications are found in nuclear power plants, aviation industry, and in automotive 
industry. 

We present a Broadband Holography Algorithm for imaging buried defects in electrically conducting 
materials as 2D-slices of the specimen. The change of conductivity is used as an indication of defects 
in homogeneous materials. The application of broadband elliptic polarized eddy current fields to 
nondestructive testing is investigated. 

A coil or a system of coils located above the test region excites a broadband polarized eddy current 
field. The depth of penetration of the electromagnetic field is highly dependent on the frequency 
of excitation, the conductivity and the permeability. The diffusion equation gives an appropriate 
description of the field distribution in metallic materials. According to the so called skin-effect, 
these fields can be regarded as highly damped ”skin-waves” which are penetrating the material 
along a limited range, and therefore testing is only possible in the near surface regions. 

We use an imaging algorithm for synthetic apertures known from microwaves [l] or acoustic waves 
[2] for the reconstruction of buried defects. We adapt the algorithm to the specific properties of 
eddy current fields. In detail, for skin-waves the phase velocity is frequency dependent (dispersion), 
phase and attenuation constants are identical. The wavelength is larger than the penetration depth 
and consequently, measurable phase shifts are very small. Furthermore, the characteristic beam 
pattern of the probe has to be taken into account because the interaction of the fields with the 
defects takes place in the near field of the probe. 

The reconstruction method can be explained by the theory of matched filtering. For this purpose 
the probe signal distribution in the receiving aperture is measured assuming that only one single 
scattering object is present (test function). Then the probe signal distribution in the aperture with 
a real defect is recorded. This result is correlated with the earlier performed test function, yielding 
a maximum value if the location where the test object has been assumed coincides with the position 
of a real scattering object. In order to reconstruct the whole image, the reconstruction algorithm 
needs the fictitious single scatterer to be moved across the whole object scene. 

Eddy current fields, excited in the medium by an arbitrary probe, are determined experimentally. 
This distribution is needed for the reconstruction algorithm. In order to achieve radial resolution 
we extend the algorithm to broadband signals (pulses, stochastic process). The concept of phase 
multiplication, proposed in [3] and [4], is applied to the Broadband Holography Algorithm for 
compensation of small phase shifts ana results in an improved radial and azimuthal resolution. 

We report on numerical experiments with simulated and measured data. The measured data is 
recorded with a test fixture in which a manipulator moves a probe in a 1-dimensional linear aperture 
over a sample containing artificial defects of different form and size. The final result shows 2D-slices 
of the tested sample. 


References 

[1] R. Karg: Multifrequency microwave holography, AEU 31(1977), pp. 150-156 

[2] H. Ermert, R. Karg: Multifrequency acoustical holography, IEEE Trans. SU-26 (1979), pp. 279-286 

[3] B.P. Hildebrand, A.J. Boland, T.J. Davis: Holographic principles applied to low frequency electromagnetic ima- 
ging in conductors, Proc. IEEE Int. Optical Computing Conf., 1983, pp. 59-66 

[4] H.D. Collins, T.J. Davis, L.J. Busse, R.P. Gribble: Eddy current phasography, Acoustical Imaging 11 (1982), pp. 
609-623 

[5] H.J. Hippier, E. Ermert, L.v. Bernus: On the Utilization of the Skin-Effect for Holographic Imaging in Nonde- 
structive Evaluation of Metallic Materials and Devices, AEU 44(1990) 4, pp. 257-265 
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Eddy Current Imaging of Defects in a Conductive Half-Space: 
Diffraction Tomographic and Generalized Solutions 
in a Wavefield Framework. 


R. de Oliveira-Bohbot, D. Lesselier and B. Duchene 
M. Nikolova and A. Mohammad-Dj afari 
Laboratoire des Signaux et Systemes (C.N.R.S.-E.S.E.) 
Plateau de Moulon, 91192 Gif- sur- Yvette Cedex, France. 


With adequate corrections and extensions, solution methods of inverse scattering problems 
can be used for Nondestructive Testing of damaged metallic structures where eddy currents are 
generated by means of external low-frequency sources. The time-harmonic diffusion equation 
satisfied by the field in — non-magnetic — metal parts is taken as a Helmholtz wave equation 
whose propagation constant is complex-valued with equal real and imaginary parts inversely 
proportional with the skin depth. Huyghens-type sources located within the defects are the 
sources of the anomalous fields observed outside the structure under testing. When successful, 
inversion of such anomalous fields provides images of the defects. But there is one severe 
challenge: the problem at hand is strongly ill-posed due to exponential decay of the probing field 
in metal (skin effect) and corresponding evanescence of fields outside. In mathematical terms, 
an integral operator with damped kernel has to be dealt with at some point of the analysis. Also, 
at least in the first step, we are restricted to linearized solutions based on the Born 
approximation (which means testing of small defects at low frequency), the inverse problem 
being nonlinear with respect to the defects 1 conductivity. 

Research on eddy-current imaging in a wavefield framework is going on since a few years 
now in our laboratory, in particular with the support of Departments of Electricite de France , 
Direction des Etudes et Recherches (EDF-DER). Present-day aspects of this research are 
outlined here with using a rather ideal two-dimensional test model: a cylindrical defect in a 
conductive half-space - a metal block in air - is illuminated by a planar source, the anomalous 
fields being observed at discrete frequencies Q) on a line parallel with the block surface. 

In the presentation the several algorithms that are currently studied are summarized and 
comparisons are carried out using results they yield from both synthetic and experimental data 
with emphasis put on their behavior vs. noise and model errors (the experiments have been led 
byN. Coutanceau at the Retour d’Experience, Mesures, Essais EDF Department): 

(i) Attenuation-matched Diffraction Tomography (DT) algorithms only account for part of 
the skin effect and provide approximate solutions but they maintain the Fourier formalism of 
the inverse scattering schemes they are issued from; DT is implemented using either data 
interpolation in the spectral space followed by fast Fourier inversion, or the method of 
Maximum Entropy Fourier synthesis, the latter being especially effective with sparse spectra. 

(ii) A DT scheme may also be employed which fully cares for skin effect though at the 
expense of an analytic continuation of the fields in the complex a) plane - this operation 
somewhat replaces an inverse Laplace transform. 

(iii) More computationally expensive inversion schemes aim at retrieving the conductivity 
in each cell of a cell-model of the defect and do not suffer from any neglection of the skin effect. 
They work in stochastic settings (e.g., Kalman filtering, simulated annealing) or in 
deterministic ones (e.g., ART, conjugate-gradient) with a priori information extensively used. 
Finally, new tools - wavelet-vaguelette decompositions (investigated with F. Brouaye) and 
Markov-random field based inversions - are introduced with first results possibly given. 
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CONSIDERATIONS ON THE DETERMINATION OF THE SPATIALLY VARYING 
CONDUCTIVITY AND PERMEABILITY OF METALS FROM EDDY-CURRENT 

DATA 

James H. Rose* 

Center for NDE 
Iowa State University 
Ames LA 50010 

ABSTRACT 

The determination of layered, one-dimensional conductivity profiles from eddy- 
current measurements has progressed rapidly in recent years. An in principal "exact" 
low-temporal-frequency method has been proposed that is based on an inverse Laplace 
transform of the impedance of an eddy-current probe as a function of the spatial 
frequency. Other, perhaps more practical methods, based on variational least-squares fits 
to the measured frequency dependence of the impedance, have also been proposed. The 
latter methods have been used by several groups to invert measurements of the 
impedance of air-core coils placed next to metals whose near-surface conductivity can be 
described by a single piece- wise constant layer (such layers might be generated by 
cladding or painting). Various groups have demonstrated the ability to determine the 
thickness and conductivity change for such single surface layer. More recently, I and 
other coworkers have successfully inverted measurements made with air-core coils placed 
near the surface of plates whose near-surface conductivity profiles vary in a smooth and 
continuous manner. We will briefly review these results for the one-dimensional case. 
Then we will discuss the prospects for two "two-dimensional" inversions. The first 
problem is to determine the magnetic permeability of a layered magnetic metal. This 
problem is two-dimensional since the permeability, }i(co, z), depends on both the depth 
and the frequency. We will also discuss the determination of the conductivity profiles of 
nonmagnetic metals that have two-dimensional spatial symmetries. 

ACKNOWLEDGMENT ; This work was supported by the Center for NDE at Iowa State 
University. 


Email jrose@cnde.iastate.edu 
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Application of Markov Models for Inverse Problems in Eddy 
Current Nondestructive Evaluation 

L. IJdpa and W. Lord 
Iowa State University 

Inverse problems in nondestructive evaluation (NDE) are involved with the reconstruction 
of defect shapes in materials, on the basis of information contained in the measured signal. 
A major difficulty in finding solutions to inverse problems in eddy current NDE lies in the 
quasi- static nature of the underlying physical phenomenon and the complex nature of the 
defect and test geometries, which necessitate the use of numerical models for solving the 
forward problem. Unfortunately these models do not provide a direct method for solving 
the inverse problem. 

This paper presents an approach that uses the numerical model for the forward problem in 
an iterative scheme for reconstructing the defects in materials. Formulating the NDE 
procedure as a mapping from the defect space to the signal space, the inverse problem 
solution is first reduced to a search in the defect space for the optimal defect profile. 
Making use of a Markov model for the defect boundary, the search space is first 
represented by a tree structure. The search space is systematically constrained using 
appropriate transition probability matrices and the solution is arrived at using a simple tree 
search procedure. 
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Parabolic Electromagnetic NDE Phenomena 

W. Lord, S. Ross, L. Udpa, and M. Lusk 
Iowa State University 


Defect characterization (imaging) capabilities of electromagnetic nondestructive 
evaluation (NDE) methods are rooted in the underlying electromagnetic field/defect interactions 
associated with the methods. Eddy current NDE phenomena are diffusive in nature and 
describable by the quasi-static or pre-Maxwellian field equations that lead to a parabolic 
governing partial differential equation (p.d.e.). This is in contrast to potential drop and flux 
leakage NDE phenomena which are static in nature and governed by an elliptic p.d.e., or 
microwave NDE phenomena which conform to Maxwell’s equations and are governed by a 
hyperbolic p.d.e. type. 

This paper discusses the imaging approaches that have been used for imaging elliptic, 
parabolic and hyperbolic NDE phenomena and introduces a new “q-transformation” between 
parabolic and hyperbolic domains, originally developed for geophysical applications, but that 
also shows promise for NDE studies. 
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AN INTRODUCTION TO LOW FREQUENCY (ULF/ELF/VLF) RADIO POLARIMETRY AND SOME 
APPLICATIONS TO THE INTERPRETATION OF NATURAL AND CULTURAL EMISSION 

SIGNATURES 

Jack Y. Dea^), Peder M. Hansen ^ \ and Wolfgang-M. Boemer^) 

(^Naval Command; Control and Ocean Surveillance Center (NCCOSC), 

Research Department, Test and Evaluation Division (NRaD), 

Communication System Development Branch (Code 832), 

Point Loma Topside, 271 Catalina Blvd, A33-1061, 

San Diego, CA, USA 92152-5000 
T/F: +[11(619)553-5781 /6457 

^University of Illinois at Chicago 
UIC-EECS Communication and Sensing Laboratory 
840 W. Taylor St., SEL-4210 
Chicago, IL, USA 60607-7018 
T/F: +[11(312)996-5480 

POLARIMETRY covers the entire electromagnetic spectrum and T.OW FREQUENCY POLARIMETRY' 
applies to the polarimetric vector (polarization) field nature of the magnetic and the electric fields within 
the ULTRA-LOW-FREQUENCY (ULF: below 3 Hz), the EXTREMELY-LOW-FREQUENCY (ELF: 3 Hz to 
3 KHz) and the VERY-LOW-FREQUENCY (VLF: 3 KHz to 30 KHz) spectral regions. The pertinent 
acoustic complementary spectral regions are identified as INFRASONIC (0.01 to 3 Hz), NEAR- 
INFRASONIC (3 to 20 Hz) and AUDIOSONIC (20 Hz to 20 KHz). These low-frequency vector 
electromagnetic and acoustic waves play an important role in coupled LITHOSPHERIC-IONOSPHERIC 
wave interactions especially in the upper ULF (above 0.01 Hz) to lower VLF (below 10 KHz) for which the 
radio waves display pronounced polarimetric effects, i.e., strong field-component dependence which may 
differ appreciably from one geographic (earth-magnetic) region to another. 

In this paper, an introduction to this existing re-emerging field is presented, putting major emphasis on 
recent findings regarding the polarimetric behavior of natural and cultural emission signature 
observations of elevated electromagnetic noise including (i) Schumann resonances generated by electric 
storm discharges (8.6, 14, 21, etc., Hz); (ii) discrete VLF signatures (5.6, 11.2, 17.2/4/6 Hz) generated 
during space vehicle passages through the ionospheric D/E/F/ layers; and (iii) signatures peaking 
strongly in the 0.1-5 Hz spectral region during lithospheric stress build-up prior to earthquake and 
volcanic distress events. The polarimetric behavior of these signatures is interpreted with the objective of 
developing inversion methods for emission source localization which then may be implemented for 
Earthquake Hazard Mitigation measures. 
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LOW FREQUENCY INVERSE SCATTERING 


Pei-Li SHEN and Shi-Ming LIN 
Northwestern Polytechnica! University 
Xian, Shaanxi 710072, P. R. China 


In order to recover the Rayleigh coefficient and the next high order nonzero 
coefficient from the backscattered signals, an inverse scattering method is pro- 
posed by Chaudhuri (IEEE Trans, on AP, Vol. AP-19,No. 3, 1981, P. 398) . 
Now, We would like to point out that there are some mistakes in the above pa- 
per. So, the values of coefficients are not true. In this paper, we will present the 
correct answer. 

Let r(t) be the returned pulse response of the object to the incident pulse Ei(t), 
and f(t) is the backscattered impulse response of the object for the given 
orientation and polarization. For the case of low frequency (ka< <1), the 
Laplace transform of f(t) can be expended as follows: 


F(s)='ZaS 


It is well known that a 0 = a,=0, and a 2 ~V, where V denotes the volume of 
scatterer. 

Now, we can prove that 

00 CO 

a 2 = c j t 3 r(t)dt / 6 j tE } (t)dt 

— oo — oo 

[ oo oo go 

— c J / 4 r{t)dt 4- 12a 2 j t 2 E i (t)dt \ / 24 J tE ( {t)dt 
— 00 — 00 J — 00 

oo oo oo 

c I t S r(t)dt + 60a 3 J / 2 (t)dt - 20a 2 J P E t (t)dt 

— oo — 00 — 00 

x £ 120 | iE t {t)dt\ 

— 00 

where c is the velocity. 


a = 

4 
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SESSION 1A2 

Monday July 12, AM 0820-1140, Baxter Lecture Hall 
ANTENNAS AND RADIATION 

Co-Chairpersons: W. A. Imbriale, JPL, California Institute of Technology, Pasadena, CA, USA, and 

J. J. Lee, Hughes Aircraft, Fullerton, CA, USA 


0820 Necessary and Sufficient Condition for the Existence of an Inflection Point 

and Its Application to Reflector Antenna Analysis 

M. H. Rahnavard, Shiraz University, Shiraz, Iran 18 

0840 Rapid Solution of Patch Resonators by a Combined Complex Image and 
Contour Integral Method 

A. A. Omar and Y. L. Chow, University of Waterloo, Waterloo, Ontario, Canada 19 

0900 Method of Linearization for the Problem of Radiation Systems Synthesis 

V. I. Dmitriev, Moscow State University, Moscow, Russia 20 

0920 Probe Compensation in Cylindrical Near-Field Scanning: A Novel Simulation Methodology 
Z. A. Hussein and Y. Rahmat-Samii, JPL, California Institute of 

Technology, Pasadena, CA, USA 21 

0940 Biquadratic Programming and Antenna Design 

K. Hirasawa, University of Tsukuba, Tsukuba, Japan 22 


1000 Coffee Break 


1020 Electromagnetic Radiation From a Series of Sharp X-Ray Pulses 

R. L. Gardner, Phillips Laboratory, Kirtland AFB, NM, USA 23 

1040 Formulation of Time-Domain Planar Near-Field Measurements 

T. B. Hansen and A. D. Yaghjian, Rome Laboratory, Hanscom AFB, MA, USA 24 


1100 Line Source Radiation in the Presence of Structures Formed by a Finite 
Number of Slotted Cylinders 

V. Veremey, Institute of Radiophysics and Electronics, Ukrainian Academy 

of Sciences, Kharkov, Ukraine 25 

1120 FDTD Analysis of a Biperiodic Phased Array of Rectangular Apertures 

J. Ren and O. P. Gandhi, University of Utah, Salt Lake City, Utah, USA 

L. R. Walker, J. Franschilla, and C. R. Boerman, Hughes Aircraft Company, 

Radar Systems Group, Los Angeles, CA, USA 


,26 
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Necessary and Sufficient Condition for the Existence of an 
Inflection Point and Its Application to Reflector Antenna 

Analysis* 


M. H. Rahnavard 
School of Engineering 
Shiraz University 
Shiraz, Iran 

In the study and analysis of surface curvature- induced shadow, which is encountered in 
dual-shaped reflectors, Glaindo- Williams surface and unstable optical resonator inflected 
point has an important role. In this paper, necessary and sufficient condition for the 
existence of inflected point is studied. The role of inflected point in analyzing an ideal 
cubic phase front is considered. 


*This research is supported by the Board of Research of Shiraz University, Shiraz, Iran. 



19 


Rapid Solution of Patch Resonators by 
a Combined Complex Image and Contour Integral Method 

A.A. Omar and Y.L. Chow 
Department of Electrical & Computer Engineering 
University of Waterloo, Waterloo, Ontario N2L 3G1, Canada 


Resonators of super-conductive microstrip patches are rapidly becoming popular because of their 
high Q -values (~1G 1 2 3 4 5 6 ). Two common analyses for such resonators are the point matching method 

[1] and the mode matching method [2]. Because the match is done over the large patch surface, 
the matrix equation resulted has frequently 180 to 400 unknowns. 

A more efficient way is to match only around an assumed magnetic wall current at the (edge) 
perimeter of the patch with 3 to 10 fold reduction in unknowns. This is the 2-B contour integral 
method of Okoshi [3] and gives reasonable results. To improve the results further, the magnetic 
wall assumption is to be relaxed to include the edge fringe fields of the patch. Martinson and Kuester 
[4] proposed a fringe correction factor depending on the substrate e r . Such factor still neglects the 
dependance on the shape of the contour of the edges. Therefore, further improvement is possible. 
The alternative for such improvement is given below. 

The resonator structure is divided into the interior and exterior regions. For the interior, we use the 
standard 2-D Okoshi formulation of [3] of magnetic current (electric voltage) on the patch perimeter. 
For the exterior, the recent 3-D complex image formulation [5] of magnetic current is used. The 
tangential electric and magnetic fields of the exterior and interior formulations are then matched 
over the patch perimeter to form a moment method matrix equation. 

In the combined method, the influence of the fringe field from the exterior is actually quite small, 
but is enough to pull Okoshi ’s solution to nearly a perfect fit with experiment. 

The studied examples are the rectangular and circular patch resonators of direct microstrip feed 
over a wide frequency range. It is evident from the formulation that the combined method can be 
applied to most patch resonators of arbitrary shape, with and without mode coupling stubs. The 
computer time is no more than 2 minutes per frequency point on a 33 MHz 80386-PC. 

[1] K.A. Michalski and D. Zheng, "Analysis of microstrip resonators of arbitrary shape," 
IEEE Trans. Microwave Theory Tech., Vol. MTT-40, No. 1, pp. 112-119, Ian. 1992. 

[2] G. DTnzeo et al., "Method of analysis and filtering properties of microwave planar 
networks," IEEE Trans. Microwave Theory Tech., Vol MTT-26, no. 7, pp. 462-471, 
July 1978. 

[3] T. Okoshi, Planar Circuits, Springer- Verlag, 1985, pp. 44-51. 

[4] T.M. Martinson and E.F. Kuester, "Accurate analysis of arbitrary shaped patch resonators 
on thin substrates," IEEE Trans. Microwave Theory Tech., Vol. MTT-36, no. 2, pp. 
324-331, Feb. 1988. 

[5] A.A. Omar and Y.L. Chow, "A solution of coplanar waveguide with air-bridges using 

complex images," IEEE Trans. Microwave Theory Tech., Vol MTT-40, no. 11, pp. 
2070-2077, Nov. 1992. 
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Method of Linearization for the Problem of Radiation Systems Synthesis 


Vladimir I. Dmitriev 
Moscow State University 
Moscow, Russia 

There is considered one multicriterial problem of the radiation system synthesis by the given 
radiation pattern. The problems of simple, optimal and quasi-optima! synthesis are investigated. 
There is developed a method of linearization for solving the problem of quasi-optimal synthesis 
for a discrete antenna lattice. It is shown how to design a computer system for discrete antenna 
systems with the given properties by means of the worked out method of the control parameters. 
Realization of such a system is described in detail and can be demonstrated on PC AT. 



